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ABSTRACT 


The class of hypersonic vehicle configurations with single stage-to-orbit (SSTO) 
capability reflect highly integrated airframe and propulsion systems. These designs are 
also known to exhibit a large degree of interaction between the airframe and engine 
dynamics. Consequently, even simplified hypersonic models are characterized by tightly 
coupled nonlinear equations of motion. In addition, hypersonic SSTO vehicles present a 
major system design challenge: The vehicle's overall mission performance is a function of 
its subsystem efficiencies, including structural, aerodynamic, propulsive, and operational. 
Further, all subsystem efficiencies are interrelated, hence, independent optimization of the 
subsystems is not likely to lead to an optimum design. Thus, it is desired to know the 
effect of various subsystem efficiencies on overall mission performance. For the purposes 
of this analysis, mission performance will be measured in terms of the payload weight 
inserted into orbit. In this report, a trajectory optimization problem is formulated for a 
generic hypersonic lifting body for a specified orbit-injection mission. A solution method 
is outlined, and results are detailed for the generic vehicle, referred to as the baseline 
model. After evaluating the performance of the baseline model, a sensitivity study is 
presented to determine the effect of various subsystem efficiencies on mission 
performance. This consists of performing a parametric analysis of the basic design 
parameters, generating a matrix of configurations, and determining the mission 
performance of each configuration. Also, the performance loss due to constraining the 
total head load experienced by the vehicle is evaluated. The key results from this analysis 
include the formulation of the sizing problem for this vehicle class using trajectory 
optimization, characteristics of the optimal trajectories, and the subsystem design 
sensitivities. 
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sonic velocity 

engine exit area 

engine nozzle area ratio 

vehicle drag coefficient 

profile drag coefficient 

vehicle lift coefficient 

specific heat at constant pressure for air 

vehicle drag 

vehicle energy height 

gravitational constant 

vehicle altitude above Earth's surface 

ratio of specific heats 

vehicle lift 

vehicle forebody length 
freestream Mach No. 
diffuser exit Mach No. 
combustor exit Mach No. 
engine exit Mach No. 
freestream pressure 
surface pressure on vehicle forebody 
freestream dynamic pressure 
vehicle heat load 

vehicle distance from Earth's center 

radius of Earth 

effective frontal area 

engine thrust 

vehicle velocity 

total vehicle weight 

effective vehicle width 

vehicle fuel flow rate 

vehicle angle of attack 

forebody nose angle 

vehicle flight path angle 

freestream density 

Earth’s rate of rotation 



I. INTRODUCTION 


A possible next-generation launch vehicle will be a fully reusable, single-stage-to- 
orbit, manned aerospacecraft. Such vehicles will most likely use liquid hydrogen fuel, be 
100 to 150 feet in length, and weigh approximately 300,000 lbs [1]. They should be 
capable of horizontal take-off and landing, accelerating to Mach 25 on orbit, as well as 
cruising at Mach 5 to 15 in the upper atmosphere. 

In order to meet the above stated mission requirements, the vehicle must utilize an 
airbreathing propulsion system for much of its trajectory, to avoid the weight penalty of 
carrying the oxidizer on board. The untested SCRAMjet engine will likely be the primary 
propulsive system. To obtain maximum propulsive efficiency, the SCRAMjet engine 
must be operated at a high dynamic pressure. However, since aerodynamic heating and 
drag also increase with dynamic pressure, the benefits of high propulsive efficiency must 
be balanced against temperature and structural constraints. Further, propulsive variables 
such as thrust and fuel flow vary with Mach number and altitude, so propulsive efficiency 
will depend on the vehicle trajectory. Similarly, both aerodynamic forces and heating are 
functions of Mach number and altitude, thus aerodynamic efficiency and structural design 
requirements (e.g., the amount of aerodynamic heating the structure can withstand) are 
trajectory-dependent as well. This demonstrates the highly interactive nature of these 
systems. As such, it should be reiterated that independent subsystem optimization will not 
lead to the optimal integrated system. 

Since these vehicles are neither conventional aircraft nor rocket, the vehicle sizing 
problem is unique. Conventional aircraft sizing, based on Breguet analysis for example, is 
not appropriate, nor is classic rocket performance analysis. Consequently, the sizing 
problem must be reformulated. A new formulation and some parametric results are 
presented in this paper. First, a baseline hypersonic vehicle configuration is defined. Also 
defined are several high-efficiency configurations, each resulting from a change in the 
efficiency of one of the vehicle's many subsystems. A minimum-fuel optimization 
problem for these configurations is then posed. After a numerical solution algorithm is 
outlined, the resulting characteristics of the baseline trajectory are detailed. Also, the 
effect of constraining the total heating load experienced by the vehicle during the mission 
will be assessed. Finally, a parametric analysis will be performed to determine the 
sensitivities between mission performance and subsystem efficiencies. 



H. DESCRIPTION OF BASELINE VEHICLE AND ATMOSPHERIC MODEL 


The vehicle geometry considered in this study is generic, but was selected to 
reflect key characteristics similar to the X-30 vehicle described in [1]. Figure 1 shows a 
simplified sketch of the baseline vehicle. The basic aerodynamic and propulsive 
modeling of the vehicle has been addressed in [4] and [5]. The configuration consists of 
a forebody/engine inlet, internal engine module, and afterbody/exhaust nozzle. 



Fig. 1. Hypersonic Vehicle Configuration. 


2. 1 Aerodynamic Characteristics 

The forebody, with forebody nose angle a n , is the primary lifting surface. If the 
vehicle angle of attack is a, under Newtonian theory the pressure on the forebody is 


p s = pao+qoo(l-cos 2(a + ctn)) 


( 2 . 1 . 1 ) 



If the pressure on the upper surface is assumed constant at p<„, the resulting lift and drag 
on the forebody are 


L= LnW e cos(a + tXi,)qoo( 1-cos 2(a + ot„)) (2.1.2) 

D = L n W e (sin(a + ctn) q«(l-cos 2(a + o^-kiooCdo ) (2.1.3) 

These equations can be arranged as 

L = qooSC L (2.1.4) 

D = qooSC D (2.1.5) 

where 

S = L n W e (2.1.6) 

C L = (l-cos(2(a + ot n )))cos (a + o^) (2.1.7) 

c d = C Do + (l-cos(2(a + a n )))sin(a + c^) (2.1.8) 


For this model, a constant Cdq of .024 has been chosen [6]. Figure 2 shows the 
aerodynamic drag polar. 



Cl 


Fig. 2. Vehicle Drag Polar. 
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The SCRAMjet engine modeling is presented in [4], [5], and [7]. The thrust and 
fuel flow characteristics of the baseline engine module are 

Thrust / (pc Ae ) = Ci F <kAk-l))- C 2 (k M J- + 1) / [M„ ffk+l)/(2(k-l))] (2.2.1) 

Fuel flow / ( p„ Ae a ) = C 3 F tfk-l)/(2(k-l)) (2.2.2) 

where 


k = 1.4 (2.2.3) 

F= 1 + ((k+l)Moo 2 )/2 (2.2.4) 

Ae = Aeoexp(-C4(Moo -1)) (2.2.5) 

Aeo = 840 ft 2 (2.2.6) 

a = 1000 ft/sec (2.2.7) 


Note that Ae is assumed to be a function of Mach number, which reflects the fact that 
the number of engine modules operating may be a Mach-dependent variable. The 
function (2.2.5) is itself a curve-fit of empirical data, with C4 taking on a value of 
approximately .2. 


For this analysis, the engine quantities have been calculated assuming fixed chosen 
values for combustor exit Mach number, nozzle area ratio, and heat input to the 
combustor. The engine dependent constants C 1 through C3 are defined below as 


= r i ih -t- m »i I *’ 

^M 3 2 + i| |l + kJ- Me 2 ) |l + M 2 2 | 



( 2 . 2 . 10 ) 



where Q and R in (3.2.9) represent heating value of the fuel and gas constant for air, 
respectively. All of the quantities in (3.2.7)-(3.2.9) are taken to be constant for a specific 
engine configuration. Therefore, the three engine constants will remain fixed over the 
trajectory. The values of the engine constants taken for this analysis are 


Ci =0.1 
C 2 = 0.03 
C 3 = 8 x 10‘ 10 /ft 


( 2 . 2 . 11 ) 

( 2 . 2 . 12 ) 

(2.2.13) 


Thrust/fPAk ) 



Fig. 3. Variation of Thrust with Mach number. 

The variation of thrust with Mach number is plotted in Figure 3. Note that this is a plot of 
Eqn (3.2.1) above multiplied by Aq/Aq 0 , thus allowing for the Mach-dependent area. 
Figure 4 shows specific fuel consumption vs. Mach number. Although thrust (lb) and fuel 
flow rate (lb/hr) are functions of Mach number and altitude. Figures 3 and 4 vary with 
Mach number only. 
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Fig. 4. Variation of S.F.C. with Mach number. 


2.3 Atmospheric Model 

For the trajectory study deiscussed later in this report, a simplified exponential 
atmosphere was used, resulting in the following model for freestream pressure poo and 
density p oo* 

(2.3.1) 

(2.3.2) 

where h is altitude and 


h ref = 150,000ft (2.3.3) 

Pref = 2-51 psf (2.3.4) 

p ref = 2.87 x 10* 6 sl/ft 3 (2.3.5) 

~ g - = 3.68 x 10-5/ft (2.3.6) 

Kgas ^ref 


P~ — Pref CXp 


— (h - h ref ) 


[Rgas Tje f 


Poo — pref exp 


(h - h ref ) 


Rgas Tref 



m. PARAMETRIC CONFIGURATION MATRIX 


Having set forth the characteristics of the baseline configuration, it can now be 
determined how to model an increase in the efficiency of the configuration's various 
subsystems in order to determine the effect of each increase on mission performance. As 
stated earlier, mission performance is measured in terms of the amount of payload inserted 
into orbit, i.e., the vehicle’s orbital mass fraction. The orbital mass fraction is the sum of 
the structural and payload mass fractions: 

W final/W initial = (W payload I W initial) + (W structure I W initial) 

So for a given orbital mass fraction, structural mass fraction trades one to one with 
payload mass fraction. Finally, the operational efficiency is maximized if the trajectory 
leads to the maximum orbital mass fraction. 

The initial and final energy levels of the vehicle's trajectory are to be specified 
(and hence known). Thus, by expressing energy as a function of mass fraction, the effects 
of the subsystem efficiencies on orbital mass fraction can be examined. Consider a 
simplified form of the energy equation as follows: 


?=£< T - D > 

«(‘ W, ) =T (w)[1'(?Il)] (assuming L = W) 



or, since dWf = -dW, 


<E - V !__l 1_ 

dWf W (s.f.c.) jL | jXj 


(3.1) 

(3.2) 

(3.3) 


(3.4) 


Equation (3.4) reveals the overall system efficiency, or the energy gain per unit of fuel 
used. In this relation one observes that the s.f.c., the thrust-to-weight ratio, and the lift-to- 



drag ratio are significant, along with the flight velocity. Obviously, the greater the system 
efficiency for a given mission, the greater the orbital mass fraction Wfmai/Winitial- 

The mission performance is therefore a function of the following subsystem 
efficiencies: 

(L/D) - Aerodynamic Efficiency 
(T/W) and s.f.c. - Propulsive Efficiencies 

In addition to the baseline vehicle configuration, three other vehicle configurations 
are generated: one representing an increase in aerodynamic efficiency, and two 

representing increases in propulsive efficiency. These configurations will be labeled as 
the high-(LZD), high-(TAV), and low-s.f.c. configurations. For each configuration, one of 
the subsystem efficiency parameters is increased while the other two are kept at their 
baseline values. Table 1 presents a configuration matrix indicating the efficiency values 
for each configuration considered. All efficiency parameters are evaluated at the nominal 
initial condition: The baseline L/D in Table 1 is defined by evaluating the lift and drag 
equations for the baseline vehicle at a = 0, while T/W and s.f.c. are evaluated at the initial 
Mach number of 2.5. Increases in each of the subsystem efficiencies are modeled as 
follows. For the high-(L/D) model, the efficiency parameters are evaluated at the same 
initial condition, but with the profile drag coefficient Cdo set t0 zero. Similarly, for the 
low-s.f.c. configuration, the engine constant C 3 for the fuel flow rate was decreased by 
10%. Finally, the high-(T/W) configuration was obtained by increasing Ci and C 2 by 
10% while also increasing C 3 by 10% to maintain the same s.f.c. as the baseline vehicle. 


Table 1 - Comparison of Baseline Configuration With 
Three Improved Subsystem Efficiency Models 


Model 


(T/W) o 

(SFC) c 

(Ib/hr/lb) 

Baseline 

2.206 

.472 

.196 

High L/D 

3.732 

.472 

.196 

HighT/W 

2.206 

.519 

.196 

Low SFC 

2.206 

.472 

.177 


As the table shows, increases in the subsystem efficiencies are as follows: 

high-(L/D): 69% increase in (L/D) 0 above baseline value 
high-(TAV): 10% increase in (T/W) 0 above baseline value 
low-s.f.c.: 10% decrease in (s.f.c.) 0 below baseline value 

As expected, each configuration shows an increase in one subsystem efficiency parameter 
and no change in the other two. 
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IV. FORMULATION OF THE TRAJECTORY OPTIMIZATION PROBLEM 

4.1 Optimization Paradigm 

The optimization problem is formulated as follows. A baseline configuration is 
selected, described in terms of what will be taken as "nominal" aerodynamic and 
propulsive subsystem efficiencies. The question is then, for this baseline system, what is 
the maximum weight, or orbital mass fraction, that can be placed in orbit? This question 
may be answered by solving for the minimum-fuel trajectory, which leads to the 
maximum orbital mass fraction and the trajectory yielding maximum operational 
efficiency for the baseline vehicle. In addition, a heating constraint on the vehicle can be 
imposed. The initial solution desired is that for unconstrained heating (i.e., the maximum 
orbital mass fraction that can be obtained without regard to heat load). Afterward, 
solutions for various specified values of maximum heat load can be obtained and 
compared to the nominal solution. The development below parallels that of [2]. 

4.2 Equations of Motion 


The single-stage-to-orbit mission is said to be accomplished if the vehicle can 
transfer from an initial energy level to an orbital energy level with non-negative final 
weight. Here energy is measured in terms of the energy height, which is the total vehicle 
energy (kinetic plus potential) per unit weight. The vehicle is treated as a point mass and 
its motion for this analysis is constrained to be planar. The equations of motion below are 
for a spherical rotating earth and are relative to a coordinate frame fixed to the earth’s 
surface [3]: 

E = r^-(Tcosa-D) + ^RsinYO)l (4.2.1) 

W b 

h = V sin 7 (4.2.2) 


7 = [( L + T sin a ) - W cos 7] + 


V 2 + (R co e ) 2 
RV 


cos 7+ 2 to e 


(4.2.3) 


W = -W f 


(4.2.4) 


The energy height, E, and the altitude above the earth center, R, are defined as 


v 2 

E = h + ^— (4.2.5) 

2g 

R = Rearth + h (4.2.6) 

The initial and final conditions are: 

h(to) = h G 
M(to) = 2.5 

Y(to)=0 

W(t 0 ) = 300,000 lb 

(4.2.7-4.2.10) 


h(tf) = hf 
M(tf) = Mf 


(4.2.11-4.2.12) 


For orbit injection, the desired values of h f and Mf will be 200,000 ft and 25, respectively. 
The heat-transfer rate near the stagnation region as given in [8] is 

0 = 865 Rn^/lWUp/p^ BTU/ft 2 /sec (4.2.13) 

where R n is the vehicle nose radius. For the baseline vehicle, the above equation becomes 

Q=Cp k 2V k l (4.2.14) 

where C is a geometry-dependent proportionality constant, and kj and k 2 are taken to be 
2.65 and 0.5, respectively. Thus the fifth equation of motion is 

Q = Q = pi/2V2.65 (4.2.15) 

with boundary conditions Q(t 0 ) = 0, and Q(tf) = Cq The parameter Cq is therefore 

proportional to the total heat load experienced by the vehicle over the trajectory. By 
selecting Cq, the total heat load may be constrained. The following values will be used 

for the known constants in the equations of motion: 


(4.2.16) 

(4.2.17) 

(4.2.18) 


Re = 4.8 x 10 7 ft 
g = 32.2 ft^/sec 

C% = 2n rad/day = 7.3 x lO^rad/sec 

The state vector for the model is then [E h y w Q]. For this model, the engine is assumed 
to be always operating at maximum thrust, so the single control is angle of attack. 

4,3 Nondimensionalization 

Before formally setting up the problem, the variables associated with the equations 
of motion (except y and a) are to be normalized with respect to reference quantities. The 
resulting nondimensional variables are: 


E = E/E ref 
h = h/h ref 
w = W/W G 

Q^Q/CPre^^ref 2 - 65 ) 

P = P/P re f 

P = P/Pref 

(4.3. 1-4.3. 6) 

where 


V = V/V ref 
R = R/R e 
T =T/(p re fAeo) 
L^Ufp^V^W,) 

D= D/( Pref^ref^LnWg) 

Wf = Wf/(C 3 p ref aAeo) 

(4.3.7-4.3.12) 


V ref = 20,000 ft/sec 

E ref =h ref +.5*(V ref )2/g 

W 0 = 300,000 lb 

h ref = 150,000ft 

Pref = 2.51 psf 

p re f = 2.87 x 10 -6 sl/ft^ 

Aeo = 840 ft 2 ' 
a = 1000 ft/sec 
C 3 = 8 x 10- 10 /ft 

(4.3.13-4.3.21) 


Thus the equations of motion in nondimensional form are: 
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E=^£ 

Eref 


X( Pre fA eo T CQS a . Pref V ref 2 LnW e p L V 
wl w 0 W 0 I 


sin y (Oe 


(4.3.22) 



^Vsiny 

“ref 


(4.3.23) 


Y= 


g 

PrefV ref2Lt)W e 

vw 

l ' W 0 


L + 


W 0 


■ T sin a - W cos y 


Rrefi 


V ref 2 V 2 + (0) e Re) 2 R 2 


RV 


+ 2 VrefOe/ 


(4.3.24) 


■yjy _ -C3 a p re f Aeo 
Wo 


exp 


)J exp[-.; 


~ h ~~ (h - 1)1 exp[-.2235(^ -1) 
L Rgas T re f I L a 


1 + 


, 2 (v^f] 

(4.3.25) 


Q= pl/2y2.65 


(4.3.26) 


The initial and final conditions are: 


h(to) — h G 
V(g = .125 
xg=o 
W(g = 1 
Q(g = o 


(4.3.27-4.3.31) 


h(t f ) = 1^ 
V(t f ) = 
Q(t f )= 


(4.3.32-4.3.34) 


where ho = 0, h f = 1.333, V f = 1.25, and Cq = CQ/(p ref 1/2 V ref 2.65). Also, E 0 = EJE k{ 
and Ef = Ef/E^. For the next several sections, all variables will be referred to without the 
bar notation. (This includes Q, which has been normalized twice, in (2.2.15) and (2.3.4).) 
Hence all quantities will be assumed to be nondimensional unless otherwise specified. 


4,4 Performance Index 
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The optimal trajectory will be that which requires the minimum fuel in order to 
ascend from an initial altitude and velocity (or energy height) to the final altitude and 
velocity in an unspecified time. Consequently, the performance index can be stated as 


min J = I Wfdt (4.4.1) 

Formulation of the problem is simplified, however, if we note that minimizing fuel is 
equivalent to maximizing final vehicle weight, W(tf). Thus the performance index can be 
restated as 

min J = -W(t f ) (4.4.2) 

This effectively changes the formulation from a Problem of Bolza to a Problem of Mayer. 
4.5 Two-Point Boundary Value Problem 


The problem as now defined is to minimize the above performance index subject 
to the above equations of motion and the specified initial and final conditions. This is a 
two-point boundary value problem, i.e., not all of the states and adjoints are known at 
either the initial time or final time. Now define the augmented performance index as 


min J = G + 


ftf rtf 

I F dt = -w(t f ) + v T ¥+ X T (f-x) dt 
j to JlQ 


(4.5.1) 


where 


¥ = [E(to)-E 0( h(t 0 )-h 0 , tfy-0, W(g.l, Q(to)-0, E(t f )-Ef, h(t f )-l*, CK^-CgF 

(4.5.2) 


A. - [Ag A h Ag Aw Aq] t 
x = [E h g W Q] T 
f = [f , f 2 f4 f 5 3 T 


(4.5.3) 

(4.5.4) 

(4.5.5) 



Note that ¥ is simply a vector of the differences between the actual and specified initial 
and final conditions placed on the states. The functions fj through f5 are the right-hand 
sides of the equations of motion in the order that they were given above. Also, define the 
Hamiltonian H = X T t The Euler-LaGrange equations— consisting of adjoint equations and 
control equations— can now be derived. The adjoint equations are 


X T = -3H/3x = -X T (3f/3x) 


(4.5.6) 


so that 


X E = -X E Of!/3E) - X h (3f 2 /3E)- X g (3f 3 /3E)- X w (3ty3E)- Xq^SE) 
X h = -XeOfj/ah) - X h (3f 2 /3h)- X g (3f3/3h)- X w (3f(/3h)- 

etc.... 


Note that the equations of motion contain terms involving quantities such as V, L, D, R, 
and T that are functions of the states. Thus, to fully express the required derivatives of f, it 
is necessary either to express these quantities in terms of the states before differentiating 
or to differentiate the quantities with respect to the states and use the chain rule. The 
control equation is 


0 = 3H/3u = 3H/3oc = X T 3f73a 

= -X E 3fj/3a - X h 3f 2 /3a- Xg3f 3 /3oc- Xy^3f^/3cc- XQ3f5/3a (4.5.7) 

Because the first two equations of motion contain L and D, which are functions of a, the 
above remark for the adjoint equations applies here. 

To derive the transversality conditions, note that there are 11 unknowns (5 states + 
5 adjoints + tf) and 8 boundary conditions. Thus, 3 new boundary conditions are expected 
from transversality. They are: 

0 = 3G/3x(tf) + (3F/3x) tf = 3G/3x(tf) + X T (tf) 

=> Xy(t f ) = 0, X w (t f ) = - 1 (4.5.8) 


0 = 3G/3t f + (F - (3F/3x)x) l tf = 3G/3t f + H l t{ 
=> H l tf = 0 


(4.5.9) 


In light of the development given in this section, the trajectory optimization 
problem can be restated as follows: solve the 5 state variables, 5 adjoints, 1 control, and tf 
that will satisfy the 5 state equations, 5 adjoint equations, control equation, and 1 1 
boundary conditions. This two-point boundary value problem may be solved, for 
example, by the "shooting" method. The method itself is described in Appendix A, while 
the application of the shooting method to this particular problem is detailed in Appendix 
B. 



V. RESULTS OF BASELINE OPTIMAL TRAJECTORY 


For the baseline vehicle, with heat load unconstrained, total flight time was 238.8 
sec., final slant range was 3.183 x 10^ ft, and final weight fraction was .5890 (i.e., 
123,300 lb of fuel expended). Time histories of the states and adjoints are given in 
Figures 5 and 6. All states are shown in the nondimensional form defined earlier except 
altitude, which is plotted in dimensional form for convenience. Note also that it is 
expressed as h - h Q , again because the orbit injection maneuver will begin at some 
nonzero altitude. 

The value of Cq for the unconstrained heating solution was found to be 
approximately 142, corresponding to a total heat soak of 5.291 x 10 5 BTU/ft 2 . Thus 
setting Cq at 142 yielded a solution for which Xq was of the order lO 7 , very near the 
necessary value of zero. 



Fig. 5. State Histories. 







lambda e 


lambda h 





Fig. 6. Adjoint Time Histories. 


Velocity and dynamic pressure are plotted in Figure 7. It is evident that before climbing, 
the vehicle must dive down below the initial altitude to a higher dynamic pressure in order 
to trade potential for kinetic energy and thus attain a more desirable flight condition. It is 
then able to ascend rapidly before levelling off slightly at the point of orbit injection. Peak 
dynamic pressure is approximately 9000 psf, which is much higher than the 1500 to 2000 
anticipated for NASP [9]. 








Fig. 7. Velocity and Dynamic Pressure Histories. 

Lift and flight path angle are plotted together in Figure 8, while drag, thrust, and 
fuel flow are plotted in Figure 9. (Note that lift, drag, and thrust are shown in the non- 
dimensional form defined earlier.) As expected, fuel flow rate is seen to increase with 
acceleration (i.e., thrust minus drag). Maximum fuel flow is extremely high 
(approximately 1600 lb/sec) compared to reasonable values for other classes of vehicles. 
Thrust turns out to be a very large quantity as well, peaking at approximately 2.1 x 10 6 lb. 



Fig. 8. Thrust, Drag, and Fuel Flow Histories. 






Fig. 9. Lift and Right Path Angle. 

To examine the effect of imposing a heating constraint, solutions were obtained for 
various values of Cq less than 140, and the resulting trajectories were investigated. 
Figures 10 and 11 show both the unconstrained trajectory and a trajectory for which the 
final heat load was constrained to approximately half that of the unconstrained 
(approximately 2.4 x 10^ BTU/ft^). Figure 10 also shows contours of constant energy E, 
dynamic pressure q, and heat-transfer rate 0 . Each of the two figures shows the 

trajectories indexed with time. Right time for the constrained case is reduced significantly 
to 84.0 sec, and final orbital mass fraction is .5864 (124,080 lb of fuel expended), a 
decrease of only about .4 % from the unconstrained result. Compared to the unconstrained 
trajectory, the trajectory with the heat-load constraint descends lower, to a greater dynamic 
pressure, and gains more kinetic energy before beginning its ascent. Note also that from 
approximately Mach 8 to Mach 17, the constrained trajectory follows a path of almost 
constant dynamic pressure at 10,000 psf. An interesting result is that the constrained 
optimal trajectory experiences much higher heating rates, but the trajectory duration is 
less than half that for the unconstrained trajectory, thus halving the total heat load. 




Comparison of Nominal ('*') & Heat-Constrained ('o') Trajectories 



lines of constant energy (ft) 

— lines of constant dynamic pressure (lb/ft A 2) 
lines of constant heating rate (BTU/sec/ft A 2) 


Fig. 10. Vehicle Trajectory-Altitude vs. Mach number. 


Comparison of Nominal ('*') & Heat-Constrained (’o') Trajectories 






VL SUBSYSTEM MISSION PERFORMANCE SENSITIVITIES 


For each vehicle configuration described in Chapter 3, the maximum orbital 
mass fraction was determined by solving numerically for the minimum-fuel trajectory. 
This was done using the shooting method (Appendix A) and following the procedure in 
Appendix B. Table 2 lists the key parametric results. Shown in this table are not only 
the maximum orbital mass fractions for each configuration, but also the mission flight 
time to achieve orbital energy, the final slant range when orbital energy is achieved, and 
the total heating load experienced by the vehicle over the entire trajectory. Also shown 
for each configuration are the same results for an optimal trajectory that is constrained 
to yield a fixed value of total heat soak. For the puposes of this analysis, a heating 
constraint of 2.48 x 10^ BTU/ft 2 was placed on all configurations. The heating 
constraint is a means of effectively decreasing the vehicle's structural efficiency. Thus 
the tabulated results below illustrate the effects of four different subsytem efficiencies 
on mission performance. 


Table 2 - Final Values of Time, Range, Total Heat Soak, 
and Total Weight Fraction for Each Model 
(For Unconstrained and Constrained Heating) 


Model 

t f (sec) 

unconst constr. 

Range (ft) 
unconst constr. 

Q f (BTU/ft A 2) 
unconst constr. 

w/w D 

unconst constr. 

Baseline 

238.8 103.9 

3.183e6 7.927e5 

5.29e5 2.48e5 

-5890 .5872 

High L/D 

222.1 98.5 

2.905e6 7.985e5 

5.10e5 2.48e5 

.5962 .5947 

High TAV 

230.3 107.4 

3.038e6 8.585e5 

4.95e5 2.48e5 

5901 .5886 

Low SFC 

234.0 102.3 

3.079e6 7.753e5 

5.30e5 2.48e5 

.6211 .6194 


The trajectories corresponding to the baseline configuration were discussed 
previously. The trajectories for the remaining cases are similar. For unconstrained 
heating, all the high-efficiency configurations achieve orbit with higher weight fraction 
than the baseline, as expected from equation (3.4). Moreover, the high-efficiency 
configurations all complete the mission faster than the baseline. As a result of the 
imposed heating constraint, the mission time for each configuration decreases 




drastically from its unconstrained value, while the penalty in orbital mass fraction is 
small— about .3% less payload for each configuration. 

Now consider the mission performance, i.e., the weight placed into orbit, 
expressed in terms of the subsystem sensitivities S 1 -S 4 by the following equation: 

A Orbital Weight = (Si) A(L/D ) 0 + (S 2 ) A(T/W) 0 + (S 3 ) A(s.f.c.) 0 + (S 4 ) AQ final 

Each sensitivity Sj is defined as the change in payload to orbit divided by the percent 
increase in the efficiency of one subsystem while holding the other subsystem 
efficiencies constant, i.e., at their baseline values. The A terms above represent the 
percent change from baseline values in (L/D) 0 , (T/W) 0 , (s.f.c.) 0 , and Qfi n ai» 
respectively. Of course, an increase in orbital weight leads directly to a corresponding 
increase in payload to orbit, if the structural mass fraction is constant and all fuel is 
expended. From the results in Table 2 and the subsystem efficiency increases presented 
in Chapter 3, assuming a gross takeoff weight of 300,000 lb, the numerical values for 
the subsystem sensitivities S 1 -S 4 for the unconstrained heating case are 

51 = 30 lb / (% increase in (L/D) 0 ) 

5 2 = 24 lb / (% increase in (T/W) 0 ) 

5 3 = 954 lb / (% decrease in (s.f.c.) 0 ) 

5 4 = 12 lb / (% increa se in Qfmal) 

If the total heat load is constrained, the subsystem efficiencies change. In this case, S 4 is 
not applicable, and the other three sensitivities become 

51 = 33 lb payload / (% increase in (L/D) 0 ) 

5 2 = 36 lb payload / (% increase in (T/W) 0 ) 

5 3 = 966 lb payload / (% decrease in (s.f.c.) 0 ) 


Note that, while the magnitudes of Si and S 3 remain about the same as for the 
unconstrained case, S 2 has increased significantly. Thus, based on this method of 


representing subystem senstivities, T/W appears to have a larger influence on mission 
performance when the heating load is constrained. The above results are summarized 
graphically in the carpet plot of Figure 12. This figure shows the variation in orbital 
mass fraction due to variations in (s.f.c.) 0 , (T/W) 0 , and (L/D) 0 . In the case of both 
constrained and unconstrained heating, the lower left comer of the plot represents the 
baseline values. 



Fig. 12. Mission Performance Sensitivities for Unconstrained ( ) and 

Constrained ( ) Heating. 


Consider again equation (3.4). If (s.f.c.) 0 , (T/W) 0 , and (L/D) 0 are known, this 
equation makes it possible to calculate approximately dE/dWf at the takeoff condition, 
or (dE/dWf) 0 . However, the values of s.f.c., T/W, and L/D— and therefore dE/dWf— 
vary widely throughout the trajectory. This suggests that it would be useful to examine 
the average values of L/D, T/W, and s.f.c. over the entire trajectory for each 



configuration. Each sensitivity above is calculated by evaluating the performance at a 
single condition (the nominal "takeoff 1 condition) and comparing it to baseline 
performance. Recall from Table 1 that for each configuration, this resulted in a higher 
value of one efficiency parameter evaluated at the "takeoff' condition, while the other 
parameters remained constant. For example, for the high-(L/D) configuration, this 
resulted in (IVD) 0 above the baseline value, and the same (T/W) 0 and (s.f.c.) 0 as the 
baseline. However, due to the high degree of coupling in the problem, (T/W) avg and 
(s.f.c.) a vg will generally not be the same as for the baseline case. This is due in part to 
the fact that decreasing the drag coefficient changes the optimal value of a. And as a 
feeds back into the equations of motion, it affects (T/W) and (L/D) as well. In other 
words, an increase in one subsystem efficiency will have some effect on other 
subsystem efficiencies over the course of the trajectory. This is borne out in the 
following tabular data. Table 3 displays (L/D) aV g, (T/W) avg , and (s.f.c.) aV g for each 
configuration, along with the average value of the parameter dE/dWf as computed by 
equation (3.4). Table 4 displays these results for the heat-constrained case. 

Table 3 - Subsystem Performance Comparison for Each Model (UNCONSTRAINED case) 


Model 


(TAV).^ 

(SFC)„ g 

(Ib/hr/Ib) 

(dE/dW f )„ g 

Baseline 

1.689 

5.108 

3.147 

.019227 

High L/D 

5.732 

5.099 

3.184 

.02090 

High T/W 

1.698 

5.128 

3.153 

.019231 

Low SFC 

1.690 

5.185 

2.852 

.02105 



Table 4 - Subsystem Performance Comparison (CONSTRAINED case) 


Model 

OVD)„ g 

07W)., g 

(SFC)„ g 

(lh/hr/lb) 

(dE/dW f)>Tg 

Baseline 

1.754 

22.943 

3.168 

.02118 

High L/D 

5.488 

21.827 

3.195 

.02139 

HighTAV 

1.749 

20.195 

3.167 

.02122 

Low SFC 

1.759 

23.821 

2.867 

.02311 


Note that both with and without the heating constraint, (L/D) avg increases significantly 
for the high-(L/D ) 0 configuration, and similarly the average specific fuel consumption 
decreases for the low-(s.f.c.) G configuration. However, the analogous situation does not 
occur for the high-(TAV ) 0 model. Indeed, (T/W) avg increases only slightly for this 
configuration in the unconstrained case, and decreases in the constrained case. This 
type of result would not be expected for a vehicle with less dynamic coupling. 

Due to the approximate nature of equation (3.4), the vehicle configurations for 
the heat-constrained case display larger average values of dE/dWf than those for the 
unconstrained case. Still, the trends exhibited by the four configurations in each case 
are clear. Consider an alternate set of subsystem sensitivites, defined as follows: 

A (dE/dW f ) avg = (Si) A(L/D ) 0 + (S 2 ) A(TAV) 0 + (S 3 ) A(s.f.c .) 0 


where 


51 = .133(% increase in (dE/dWflavg) / (% increase in (L/D) 0 ) 

5 2 = .0208(% decrease in (dE/dWf^yg) / (% increase in (TAV) 0 ) 

5 3 = .948(% decrease in (dE/dWf)^ vg ) / (% decrease in (s.f.c.) 0 ) 


for no heating constraint, and 




51 = .0144(% increase in (dE/dWf)avg) / (% increase in (L/D) 0 ) 

5 2 = .0189(% decrease in (dE/dWf)avg) / (% increase in (T/W) 0 ) 

5 3 = .91 1(% decrease in (dE/dWf^ V g) / (% decrease in (s.f.c.) 0 ) 

with heating constraint. Again, due to simplifications in the derivation of equation 
(3.4), the computed values of (dE/dWf) aV g in Tables 3 and 4 are approximate. 
However, the general trends for the Sj's are in line with those of the Si's defined 
previously—specifically, the sensitivity S 3 is much larger than Sj or S 2 both with and 
without heating constraint. Thus the data in Tables 3 and 4 offers some insight as to the 
effect that a change in one subsystem efficiency has on the overall system efficiency, as 
well as on the other subsystem performance parameters, over the entire trajectory. Both 
sets of sensitivities defined in this section suggest that, from a vehicle designer's 
perspective, vehicle fuel efficiency is the most critical design aspect with regard to its 
effect on overall system efficiency. As stated in Section 3, maximum overall system 
efficiency translates directly into maximum mission performance, i.e., maximum 
payload to orbit. 



vn. CONCLUSIONS 


Single-stage-to-orbit (SSTO) vehicles have unique characteristics. Like an aircraft, 
they rely on significant aerodynamic lift and have smaller thrust-to-weight ratios than 
rockets. But like rockets, their mission is to maximize payload into orbit, and propulsive 
and structural efficiencies are extremely important. The orbit-insertion mission here is 
more one of climbing in the atmosphere than cruising like an aircraft or a ballistic launch. 
But the sizing of SSTO vehicles leads to a unique problem. The sizing problem has been 
formulated for a generic hypersonic vehicle model utilizing an optimization paradigm, 
built around minimum-fuel optimal trajectories, and defining mission performance in 
terms of the maximum amount of payload placed into orbit. The optimal trajectory 
involves a zoom/climb, reminiscent of classical aircraft climb performance. Although 
orbital mass fraction for the mission is approximately 60% for all the configurations 
considered here, results indicate that this value is a strong function of L/D and s.f.c., and 
depends on the T/W ratio as well. In addition, a heating model has been incorporated into 
the problem as a means of constraining the heat load experienced by the vehicle. 
Although a heating constraint clearly can reduce the maximum orbital mass fraction for a 
given vehicle, trajectory shaping via optimization can significantly reduce this deleterious 
effect. 
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A. NUMERICAL SOLUTION METHOD 


A.1 Shooting Method 

The shooting method for solving two-point boundary value problems involves 
guessing the unspecified initial conditions, integrating forward, and updating the guesses 
(via a Newton-Raphson routine) based on the error between the actual and desired 
terminal conditions. The method as applied to this particular problem proceeds as 
follows: 

Define the states and adjoints together as X = [Xj XJ 7 , where Xj = [E h y W Q] T are the 
elements of X that are specified at the initial time and X 2 = [Xjr Xy X w Xq] t are the 
elements of X that are unspecified at the initial time. (Note that for this problem X] is 
comprised entirely of states and X 2 entirely of adjoints, although generally this may not be 
the case.) Thus the state and adjoint equations can be written as 

X = F, where F = [f -3H/a X ]T (A. 1.1) 

Define M(t) = aX(t)/dX 2 (t 0 ). From this we see that 

M(t 0 ) = ax(t 0 )/ax 2 (t 0 ) = [0 i 5 ]T ( a.i.2) 


and 


m = ax (t)/ax 2 (t 0 ) = aF/ax 2 (t 0 ) = aF/axa)*ax(t)/ax 2 (t 0 ) = aF/ax(o*M 

(A. 1.3) 

With the initial condition on M (4.1.2), Equation (4.1.3) can be integrated together with 
(4.1.1) at each time step. Now define 

¥ m (t) = [E(t)-E f , h(t)-h f , Xy(i)-Xy ( , X w (t)-X wf , H(t)-H f f (A. 1.4) 

where Ef, hf, etc., are the final conditions specified for the states, adjoints, and 
Hamiltonian. So 


¥ m (t f ) = [E(t f )-E f , h(t f )-lf, ty(t f U W (tf)+l,H(tf)]T 


(A. 1.5) 


Note that 

a'F m (t f )/ax 2 (t 0 ) = 34> m (t f )/ax(t f )*ax(t f )/ax 2 (g = a'f m (t f )/ax(t f )*M(t f ) 

(A. 1.6) 

After (A.1.1) and (A. 1.3) are integrated up to tf, M(tf) is known and a4 / m (tf)/3X(t f ) can be 
calculated. Thus 3'F m (tf)/ax 2 (t 0 ), a 5x5 matrix, is known. 

To solve for X 2 (t 0 ) (which for this problem consists of the initial adjoints X(g) 
that will drive 'I'mCtf) to its desired value of [0], an initial guess for X 2 (g is made and the 
state equations for X and M are integrated from ^ to tf. The Taylor series expansion of 
'F m(tf) as a function of X 2 (g gives 

'F m ( t f)i+1 = 'Fm(tf), + O i F m (t f )/aX 2 (t 0 )) i * AX 2 (to)i (A. 1 . 7 ) 

where the desired value of v F m (t f ) i+1 is [0], 4^(1^ and (a v F m (t f )/aX 2 (g) i are computed 
from the integration, and AX 2 (t 0 ), gives the amount by which to increment X 2 (to) in order 
to reach the desired M^tf). So the initial guess for X 2 (g is updated according to 

X 2( t o)i+l = X 2( t o)l + AX 2( l o)i 


= X 2 M + O v F m (t f )/ax 2 (t o ))-l 0 F m (t f ) desired - 4' m (tf)i) 


= X 2 (to)i - (a^ m (t f )/ax 2 (t 0 ))-^ m (t f ) i (A. 1.8) 

Because changes in 4' m (tf) can be very sensitive to changes in X 2 (t 0 ), this equation is 
modified by 


X 2 (to)i + l = XzCg - o(d'F m (t f )/ax 2 (t 0 ))-l'F m (t f ) i (A. 1.9) 

where a is a step size parameter between 0 and 1. Use the updated X 2 (to) to integrate 
again from t G to tf. Iterating in this manner will cause X 2 (t 0 ) to converge to X 2 *(g 
corresponding to the optimal solution. There are a number of possible critieria to 
determine convergence, such as ll v F m (tf) i ll < e or IIX 2 (t 0 ^ +1 - X^t^ll < £. 
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Two issues still to be addressed are determination of tf and determination of the 
optimal control 0 ^. 

A.2 Determination of t f 

Six boundary conditions have been specified or derived at the final time. Five of 
them are used in the definition of 4^(1) in order to determine the five initial adjoints 
X 2 (t 0 ) via Equation (A.1.9). The sixth, Q(t f ) = Cq, is used as a stopping condition (i.e., 
the terminal time tf has been reached when the value of Q reaches Cq). Any of the final 
conditions can be used for this purpose. However, as can be seen from the fifth EOM, Q 
always increases monotonically from zero. This guarantees that the value of Q will reach 
Cq at some point and will reach it only once, making it a good candidate for a stopping 
condition. 

This raises the question of how to choose Cq. Different minimum-fuel trajectories 
can be achieved for various choices of Cq, each representing a certain degree of heating 
constraint. Initially, it is desired to find the minimum-fuel trajectory with no heating 
constraint. From the derivation of the Euler-LaGrange adjoint equations, it is readily seen 
that Xq = 0. That is, because f is not a function of Q, the quantity }J(dfJdQ) is zero. 
Therefore, Xq is constant throughout the trajectory. Since Xq indicates the sensitivity of J 
to changes in Q, the "unconstrained heating solution" will be the minimum-fuel trajectory 
corresponding to Xq = 0. Of course, the value of Cq that will yield this solution is not 
known a priori. Yet, as long as a solution can be obtained for some arbitrary value of Cq, 
it is likely that the unconstrained heating solution can be found by repeatedly perturbing 
Cq and recomputing a solution until Xq is near zero to within some tolerance. Once the 
unconstrained solution is obtained, Cq can then be varied to yield heat-constrained 
solutions. 

A.3 Determination of oc opt 

The optimal control, which satisfies the control equation 9H/9a=0 (2.5.9), is 
determined at each time step by the Newton iteration formula 

a i+1 = cq - [0H/3a)/02H/aa2)] a=ai (A.3.1) 

The initial guess for a opt isl° at t^ and thereafter is the previous value of a opt . There 
are a number of possible critieria to determine convergence, such as l0H/5odl < e or 
H04+1 - cqll < £• 


B. SOLUTION PROCEDURE 


B.l Baseline Configuration 

Due to the tight coupling and nonlinearities inherent in the state equations, solution 
of the full mission trajectory by the above method tends to require extremely accurate 
guesses of the initial adjoints A,(t 0 ). In addition, because the adjoints are contrived 
parameters and not physical quantities, it is difficult to determine what the magnitude of 
these parameters should be. The problem was solved by means of an extrapolation or 
"bootstrapping" procedure, which is described as follows: 

1) Pick an "intermediate" trajectory, i.e., one that begins at the initial flight 
condition given in Chapter 2, but terminates at an altitude and velocity much lower than 
orbital values. An estimate of the value of Cq (the specified terminal condition on Q) for 
this trajectory is obtained by guessing X(to), integrating forward from the initial state until 
one of the specified terminal values (such as altitude or velocity) is reached, and noting the 
final value of Q. For several different guesses of A,^), the values of Q(tf)-and the other 
final state values-give an indication as to the size of Cq required. The numerical solution 
method of the previous section is then applied to the intermediate trajectory. (Note that in 
picking an intermediate trajectory, the final conditions on E, h, and Q are different from 
their values for the full trajectory. However, by transversality, the conditions on A.y(tf), 
A. w (t f ), and H(tf) remain the same.) 

2) Once X,(t G ) corresponding to the optimal solution (call this X*(t 0 ) ) is found for 
this short trajectory, the next step is to obtain a solution for a longer trajectory (i.e., a 
higher energy level) by extrapolation. One method of extrapolating from one trajectory to 
the next is to begin with the current X*(t G ) and integrate past the terminal conditions of the 
current trajectory to a higher energy level. The values of E, h, and Cq at this energy level 
become the new E f , h f , and Cq, and the current X*(t D ) becomes the initial guess of X.(to) 
for the new trajectory. The numerical solution method is then applied to the new 
trajectory. 

3) As this process continues, the specified terminal conditions for each intermediate 
trajectory are made to approach those for the full (desired) trajectory until a solution for 
the full orbit injection mission trajectory is reached. 


As the above solution procedure was followed, it was found that for each 
intermediate trajectory, the smaller the change in energy level from the previous 
trajectory, the faster the solution method converged. It should also be noted that the 
minimum-fuel solution for each intermediate trajectory need not be the one corresponding 
to unconstrained heating load for that trajectory. At each intermediate step, finding a 
minimum-fuel solution for some value of Cq-not necessarily the one for unconstrained 
heating-made it possible to march toward the final trajectory by extrapolation. Only in 
the case of the final trajectory is the particular value of Cq of importance. 

B.2 High Efficiency Configurations 

In solving the minimum fuel trajectory for each high efficiency configuration, the 
full trajectory solution for the baseline model was used as a "starting" solution. The 
general procedure is here described for the low-s.f.c. configuration: The engine constant 
C 3 was first decreased by approximately 1%. Using A.*^) for the baseline trajectory as 
the initial guess for A(t 0 ), a convergent solution was obtained for the perturbed 
configuration by the shooting method. Then, using the current values of A*(t 0 ) as an 
initial guess for A(t 0 ), a solution was obtained for C 3 decreased slightly further. This 
process continued, perturbing C 3 a little at a time and re-optimizing, until the constant was 
10% below its baseline value. Similarly, in the case of the high T/W configuration, the 
engine constants Cl and C 2 were increased slowly from their baseline values until they 
were 10% higher; and for the high L/D model, the Cdo was decreased slowly from its 
baseline value to zero. 


